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Factorization of numbers with the help of Gauss sums relies on an intimate relationship 
between the maxima of these functions and the factors. Indeed, when we restrict ourselves to 
integer arguments of the Gauss sum we profit from a one-to-one relationship. As a result the 
identification of factors by the maxima is unique. However, for non-integer arguments such 
as rational numbers this powerful instrument to find factors breaks down. We develop new 
strategies for factoring numbers using Gauss sums at rational arguments. This approach may 
find application in a recent suggestion to factor numbers using an light interferometer [V. 
Tamma et al., J. Mod. Opt. in this volume] discussed in this issue. 
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1. Introduction 

Newtonian mechanics originates only from those trajectories in spacetime, which 
make the action an extremum. When we free ourselves from this restriction and al- 
low all trajectories, we arrive at quantum mechanics provided we introduce another 
selection criterion: although all trajectories contribute in a democratic way, that is 
with equal weight, each trajectory carries a phase determined by the action along 
this path. In this way we obtain the path-integral formulation (Q, [j]) of quantum 
mechanics pioneered by Richard P. Feynman Qj. More freedom combined with a 
new rule in this enlarged realm gives birth to a new theory, which transcends the 
original one. 

On a much more elementary level, the factorization of numbers with the help 
of Gauss sums(0, Q) displays a similar phenomenon. So far, we have restricted 
our analysis of these sums to integer arguments. This requirement has provided 
us with a tool for the unique identification of the factors. In the present paper 
we allow more freedom and consider these sums at rational numbers. Needless to 
say, we immediately loose the possibility of identifying factors. However, when we 
introduce new rules for extracting factors the enlarged space of rational numbers 
enables us to factor numbers more efficiently. 



* Corresponding author. Email: sabine.woelk@uni-ulm.de 

1 It is interesting to note that already in 1924, that is before the advent of modern quantum mechanics 
introduced by Heisenberg and Schrodinger, Gregor Wentzel proposed the path integral in a paper entitled 
"Zur Quantenoptik". Unfortunately, he never followed up on his idea. 
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1.1. Formulation of the problem 



Factorization of large integers is an important problem in network and security 
systems. The most promising development in this domain is the Shor algorithm(0), 
which uses the entanglement of quantum systems to obtain an exponential speed 
up. 

Recently a different route towards factorization was proposed. Based solely on 
interference, an analogue computer evaluates a Gauss sum. Several experimental 
implementations have been suggested (0, Q, 0)- In the meantime Gauss sum 
factorization has been demonstrated experimentally in various systems ranging from 
NMR techniques (f^-[7^). cold atoms (fiil-173). ultra short laser pulses (1^, \ldi ) to 



classical light interferometry (fTj , 11). The largest number factored so far had 17 
digits. 

All experiments performed so far implement the truncated Gauss sum 



m=0 ^ ' 



where iV is the number to be factored and the integer £ represents a trial factor. 

The general idea of this approach is to find an observable of a physical system 
such as a spin, or the internal states of an atom, which is given by such a Gauss sum. 
Obviously a requirement for the successful application of Gauss sum factorization is 
that the parameters £, M and N can be chosen at will. In this way, we implement 
an analogue computer, that is the system calculates for us the truncated Gauss 
sum A^\t). Unfortunately the experiment does not provide us directly with the 
factors but only with a "yes" or "no" answer to the question if £ is a factor. Indeed, 
the result of the experiment is: 

"£ is a factor of N" or "£ is not a factor of N" 



As a consequence, we have to check every prime number up to y/N. Most of the 
experiments performed so far in the realm of Gauss sum factorization had to check 
every trial factor £ individually. However, there already exists an experiment (73, 73) 



with classical light, where the Gauss sum A^\tj is estimated simultaneously for 
every trial factor £. Here the number £ is encoded in the wavelength of the light 
and we deduce the factor of N from the spectrum of the light. In this way, we have 
the opportunity to obtain all factors simultaneously. 

However, as often in life, every opportunity is acc omp anied by new challenges^. 
Indeed, the spectrum in the experiment of Ref. (fTi . 11 ) is not represented by the 



Gauss sum at integer values of £ only. It is determined by the Gauss sum evaluated 
along the real axis. So far, we have only used the Gauss sum of integers £ to factor N. 
As a consequence we throw away a lot of information contained in the full spectrum. 
The goal of the present paper is to make the first step towards factorization with 
Gauss sums at real numbers: we discuss factorization with truncated Gauss sum at 
rational arguments. 



1 One of the referees has speculated in his report that if each challenge would carry a phase we may find a 
new quantum mechanics of life. Indeed, Pascual Jordan ll 1 3i ) and Erwin Schrodinger i 1 <A ) have argued that 
life might be intimately connected to quantum mechanics. 



October 25, 2012 0:33 Journal of Modern Optics Gauss_Sum 



3 

1.2. Outline 

Our paper is organized as follows. In section [2] we briefly summarize the theoretical 
and experimental work on Gauss sums. Then we analyse in section [3] the behavior 
of the truncated Gauss sum A^f for rational arguments. Our main interest is to 
obtain an answer to the following question: is it possible to extract information 
about the factors of N from peaks of the Gauss sum located at rational arguments. 
We show, that the answer is an emphatical "yes!" by providing two methods. We 
conclude in section [4] by a brief summary of our results. 



2. A family of Gauss sums 

Gauss sums come in a large variety of forms and have already been analyzed (0) 
in their possibilities to factor numbers. In this section, we give a brief introduction 
into the different types of Gauss sums and how to obtain factors from them. We 
conclude by summarizing the main theoretical and experimental results obtained 
so far. 



2.1. Continuous Gauss sum 

The Gauss sum 

oo 
m=— oo 

with the weight factors w m depends on the continuous argument £ and contains an 
infinite number of terms. We obtain the factors of N by testing if a maximum of this 
function is located at an integer value of £ as shown in FigJT] The sum emerges (jffZl ) 
in the context of wave packet dynamics in the form of the autocorrelation function, 
or the excitation of a multi-level atom with chirped pulses tt2fi ). 



2iri I m + 



m 
~N 



(2) 



I s 91 (01 




Figure 1. Factorization of TV = 91 = 7 ■ 13 using the continuous Gauss sum Sjv(£) defined by Eq. J2j- On 
the top we present |Sgi(£)| 2 in its dependence on £ for 1 < £ < 16. On the bottom we magnify |S9i(§)| 2 
in the vicinity of candidate primes. The pronounced maxima at the factors 7 and 13 are clearly visible. In 
contrast, at non-factors the signal does not show any peculiarities. Here, we have used Gaussian weight 
factors Wm centered at m = of width 5m = 10. 



October 25, 2012 



0:33 



Journal of Modern Optics Gauss _ Sum 



2.2. Gauss sums with integer arguments 

The Gauss sum 



Sn(£) = ^ w m ex.p 



m=—oo 



2 * 

2in m — 

N 



(3) 



follows from the continuous Gauss sum Eq. (|2|), when we restrict the argument £ to 
integer values t. Since exp [2irim£] = 1 the term linear in m disappears and only 
the fraction ijN survives. 

In Ref.(0) we transform Sn(£) defined by Eq.(|3j) into the standard Gauss sum 



N-l 



m=0 



2 l 

2iti m — 
N 



(4) 



In contrast to the Gauss sums SV(£) and Sn{£) we now have a finite sum, which 
contains N terms. This simplification emerges^ from the fact, that the phase 
factors 7 m = exp \2-Kim 2 £ / N 1 ^ are periodic in m with period N, that is ^ m +N = 7m- 
For the factorization of numbers, it is useful to introduce the scaled square of the 
standard Gauss sum 



1 



g N {e) = -\G{l,N)Y 



(5) 



The standard Gauss sum G(£,N) can be calculated analytically (22). Due to the 
division of G(£, N) by N we find the result 



g N (£)=gcd(£,N) (6) 

which gives us three possibilities to find a factor q of N: 

(i) If £ = n ■ q, that is £ is an integer multiple of a factor, <?7v(n • q) is larger than 
1, whereas for all other £ the function Qn{£) is equal to unity. 

(ii) The value of Qn{£) is a factor of N, if I is an integer multiple of this factor, 
that is gjv(n ■ q) = q. 

(iii) The function Qn{£) is periodic with the period q. 

In Figj2]we illustrate this technique for N = 91 = 7 • 13. We note that the function 
<7gi (£) displays maxima at integer multiples of 7 and 13 with the values 7 and 13. 

2.3. Reciprocate Gauss sum 

The Gauss sum 



A 



!v M) < 



1 M 

m=0 



exp 



2 N 
2irim — 



(7) 



which is most popular in the context of experiments is again defined for integer 
arguments but uses the fraction N/ £, that is the reciprocate of the ratio appearing 
in G(£,N). Morover, in the summation extends only over M terms, which is 



October 25, 2012 0:33 Journal of Modern Optics Gauss_Sum 



13 
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7 • • • • • 



7 14 21 28 35 

Figure 2. Factorization of N = 91 = 7 ■ 13 using the scaled square of the standard Gauss sum gjv(^) = 
\G(e,N)\ 2 /N denned by Eqs.gJ and (BJ. In contrast to FigQ] were we had the continuous variable £ we 
now restrict ourselves to discrete arguments I. We recognize dominant maxima at integer multiples of the 
factors 7 and 13. Moreover, we note the relation gpf(n ■ q) = q for integer n and the factor q of N. Indeed, 
the value of t/jv a t an integer multiple of a factor is the factor. 

1 h • • • 

1/V2 - 

. . .* . »• .*•.*»•• .. •» * % .*«*. .« .. ^ 

7 13 91 

Figure 3. Factorization of N = 91 = 7 ■ 13 using the truncated Gauss sum A^f\ti defined by Eq.Q. 
This function assumes the value unity only for factors. For non-factors, the function falls below 1/V2 
provided M is large enough. Here we have used M = 8 for the evaluation of the Gauss sum. 

not necessarily identical to N. For this reason, this Gauss sum carries the name 
"Truncated Gauss sum". 

The truncated Gauss sum A^\i) is of special interest, because it needs only a 
few summands to distinguish between factors and nonfactors. Provided the min- 
imum number of terms grows (0, 0) with </N, the absolute value of A { ^f\t) 
is smaller than l/\/2 for nonfactors and equal to unity for factors, as illustrated 
in Fig. [3] For this reason all experiments on factorization with Gauss sums were 
performed with 

2.3.1. Experiments 

Many experiments are based on the interaction of laser pulses with two-level 
systems. Here the fraction N/l is encoded in the phases of the pulses. After the 
sequence of pulses the occupation probability of the excited level is given by the 



truncated Gauss sum. Such experiments were performed with cold atoms () 1 if ) and 
liquids using Nuclear Magnetic Resonance (f^-U^). 

Other experiments relied on a Bose-Einstein-Condensate in an optical lattice 



() 12 , 1 1S \. From the measured momentum distribution, the average energy, which 
is given by the truncated Gauss sum, can be obtained. This technique is very 
interesting, since it was claimed that the probability for higher momentum itself, 
can be a more sensitive measurement. 

All experiments so far share the problem (0), that they have to calculate the 
ratio N/l to determine the phases of the pulses. Fortunately there already exist 
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some proposals for experiments with light, which do not suffer from this problem 



some proposa 

QBQ. The first proposal \U) uses a sequence of interfering femto-second 
laser pulses with identical phases but variable separation between the pulses. The 
recent experiment ([73 , 11 ) employs an optical interferometer with variable optical 



path lengths. In both of the proposals, the truncated Gauss sum is calculated not 
only for integer values, but also for non-integers. This aspect has motivated us, to 
investigate the truncated Gauss sum at rational arguments as discussed in the next 
section. 

2.3.2. Methods to avoid ghost factors 

The truncated Gauss sums have the problem, that for some non-factors, their 
values are close to unity provided we do not take into account enough terms in the 
sum. We have called these factors ghost factors. 

There exist two techniques to improve factorization with truncated Gauss sums 
and in particular, to decrease the number M : (i) the first one relies on the expo- 
nential sums 



M + 



1 M 



m=0 



2ni m? — 



(8) 



Indeed, we avoid (73, 2£) ghost factors, if M is larger than, or of the order of ^/N 
(ii)The second method ([73, 

0) 

uses M randomly chosen phases instead of M 
consecutive phases. With this method the minimum number M of terms necessary 
to factor a number grows only logarithmically. In this way it was possible to factor 
a 17-digits number with only 10 pulses. 



3. More freedom:Truncated Gauss sum at rational arguments 

So far, we have only considered the truncated Gauss sum at integer arguments 

£. This restriction was motivated by the fact, that A.^ {£) is equal to unity, if and 
only if the integer trial factor £ is a factor of N. In the present section we show 
that we can also gain information about the factors p^ of N from Gauss sums at 
rational numbers. 



3.1. Three classes of rational arguments 

Such an extension of the theory has been made necessary by a recent experiment 
(11 a . Lo) where all trial factors are checked simultaneously. This experiment involves 
the truncated Gauss sum not only for integer but also for rational numbers. 
Indeed, the generalized truncated Gauss sum 



< ) (0- F VlEe X p(2^m 2 f) (9) 

m=0 ^ ? ' 

for the continuous variable £ assumes the value A^f^ (£) = 1 for rational numbers of 
the form 1/s, N/s and Pk/s with s integer and p/% factor of N. The resulting peaks 
have the information about the factors encoded either in their location, or in their 
frequency of appearance, or do not contain any information at all. Obviously the 
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maxima at £ = Pk/s contain information since their locations are proportional to 
the factors p^- Since the peaks at £ = 1/s are independent of N, they cannot give 
any hints about factors of N. Unfortunately, the analysis of the peaks at £ = N/s 
is not that straight forward. Nevertheless, we show in the next section, that it is 
possible to obtain information about the factors 



3.2. New selection criteria 

The problem with the truncated Gauss sum at rational arguments as a tool to 
factor numbers originates from the existence of additional peaks located at £ = 1/s 
and £ = N/s. On first sight, these maxima do not seem to give us any information 
about the factors of N. However, in the present section, we outline two methods, 
that allow us to extract the factors from these additional peaks. 

3.2.1. Reduction and discretization of search interval 

If the Gauss sum would not contain the peaks at £ = 1/s or £ = N/s, we could 
determine the factors immediately from the maxima at p/s. All we have to do in 

this case, is to search for a datapoint, where A.^ (£) = 1. The value £ = Pk/s will 
then give us a factor pk of N . 

Hence, we have to develop a strategy how to eliminate unwanted peaks and/or 
identify those peaks which contain useful information. 

From the peaks at £ = 1/s we cannot learn anything, since 1/s does not contain 
any information about the factors. Fortunately, it is easy to ignore those peaks, 
since they are all located in the domain £ < 1. For this reason, we restrict our 
domain of interest to arguments 1 < £. 

We now have to select from the forrest of peaks arising for 1 < £ the ones, 
which carry information about the factors. For this purpose, we discretize the curve 
obtained from the Gauss sum A.^ with continuous argument £ by considering 
values of £, which are integer multiples of the minimal step size 1/sq. Moreover, we 
confine the measurement range to 1 < £ < y/~N /sq. This domain must contain peaks 
at £ = Pk/s, but there will be no peaks at £ = N/s. The following consideration 
explains this observation. 

We assume, that at a given multiple integer I of the step size 1/sq, there exists 
a peak, which can be attributed to the ration N/s, where N and s do not share a 
commen factor. As a result, we have the identity 

or i=^N (10) 
S SQ s 

Since s does not share a factor with N and I is an integer, s has to be a divisor of sq. 
However, this fact implies the inequality s < so and from Eq, (|10p we find N < I. 
This inequality is in contradiction with our assumption of the domain 1 < £ < 
which provides us with the inequality £ = £/sq < y/N /sq, that is £ < y/N. 

The reduction of the interval together with the discretization of the arguments 
of the truncated Gauss sum ensure that all peaks appearing in the so-obtained 
diagram carry the information about the factors of N. Figure [4] illustrates this 
method. 

We conclude this discussion by recalling that in the standard approach of factoring 
numbers using Gauss sums at integer values, we have to search for factors at prime 
numbers smaller than y/N. In the worst case, we have to test all prime numbers up 
to y/N. In the present approach, we search for factors at integer multiples of 1/sq 
in an interval up to y/N / sq. As a result, the present method scales in the same way 
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1/V2 



1 13 I 2 13 3 

10 5 5 



Figure 4. Factorization using the truncated Gauss sum Ajj (£) at rational arguments § = r/s, with r 
and s integer, illustrated by the example N = 91 = 7 • 13. Here we have confined ourselves to the range 
1 < £ < 3 and have chosen a step size 1/so = 1/10. Despite of the fact that there are no factors in this 
domain, the Gauss sum assumes the value unity at the positions £ = 13/10,7/5 and 13/5, corresponding 
to the maxima at § = Pk/s where p^, is a factor of N. This analysis immediately suggests the factors 7 
and 13. Due to the reduced domain of interest and the discreteness of the observation points, we do not 
observe the peaks from arguments £ = 1/s or £ = N/s except the one located at £ = 1. 



as the standard one. However, the interval in which we perform the search is much 
smaller. It is in this sense that factorization with Gauss sums at rationa arguments 
might be more efficient than the corresponding one at integers. Moreover, it is 
appealing to use the full information contained in the Gauss sum. 

3.2.2. Degeneracy of ratios 

We now propose a second method to take advantage of the additional peaks in 
the Gauss sum at rational numbers. For this purpose we analyse the series N/s 
and pk/s and note, that pk/s is contained in N/s. In other words, because pk is 
factor of N, some numbers N/s can be expressed by Pk/s'. The number of different 
representations of the fraction N/s depends of course on the factors of N. 

In figure |5] we show for the example ofiV = 3- 7 = 21 the degree D of degeneracy 
of these different representations of the same number, providing information about 
the factors. Indeed, the periods contained in D are the factors. 

We conclude by briefly addressing the question of how to implement the degree 
of degeneracy of ratios in a quantum system. Many ideas offer themselves. In this 
context it suffices to name at least one. For this purpose we recall the phenomenon 
of quantum carpets Here structures appear in the spacetime representation 

of a Schrodinger particle moving between two hard walls. Indeed, the design of a 
quantum carpet has its origin in the degeneracy of the so-called intermode traces 



§2% . It is only a small additional step to conjecture that the degeneracy of the 
ratios manifests itself in different steepnesses of the canals. However, a more de- 
tailed analysis goes beyond the scope of the present paper and we refer to a future 
publication. 



4. Summary 

The idea of factoring numbers using truncated Gauss sums relies on the fact, that 
these sums assume the value unity if and only if the test factor is indeed a factor. 
However, this one-to-one relation is only true when we restrict the arguments of the 
truncated Gauss sum to integer values. When we give up this restriction, we loose 
the possibility of uniquely identifying factors. In this case the truncated Gauss sum 
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Figure 5. Factorization ofiV = 21 = 3- 7 with the help of the number D of degeneracies of the truncated 
Gauss sum A^ 1 ^ at rational arguments. For N = p ■ q this function assumes the value unity at the 
rational numbers p/s,q/s and N/s, marked on the corresponding axes by rhombs, squares and full dots, 
respectively. Fractions, which are connected with factors allow several representations. Indeed, the ratios 
21/6 and 7/2 display the degree D = 2 of degeneracy as indicated in the top by the full dot in the square. 
Only the ratio 21/21 is identical to 3/3 and 7/7, giving rise to D = 3. The degree of degeneracies displays 
a double periodicity determined by the two factors. For example, the squares with the dots have the period 
3, whereas the rhombs with the dots show the period 7. 



takes on the value of unity at many non-integer values, which obviously cannot be 
factors. However, this additional wealth of maxima also opens up new possibilities 
of factoring numbers. In the present paper we have introduced two methods, that 
allow us to obtain with the help of these on first sight useless peaks additional 
information about the factors. 
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